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Abstract 

The observation of three events for the decay pfi^fi^ with a dimuon invariant mass of 

214.3 lb O.SMeV by the HyperCP collaboration imply that a new particle X may be needed to 
explain the observed dimuon invariant mass distribution. We show that there are regions in the 
SUSY-FCNC parameter space where the in the NMSSM can be used to explain the HyperCP 
events without contradicting all the existing constraints from the measurements of the kaon decays, 
and the constraints from the — mixing are automatically satisfied once the constraints from 
kaon decays are satisfied. 
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I. INTRODUCTION 



Recently, there are a great deal of interest [l|, l2|, y, 14, l5|, l6|, 17|, l8|, l9|, Il0| in the interpretation 
of the observed three events for the decay S+ —>■ pfi^fi'" with a dimuon invariant mass of 
214.3 ± O.SMeV, which were reported by the HyperCP Collaboration The branching 
ratio based on the three events for this process is (8.6^^:^ ± 5.5) X It has been 

argued that in the framework of the Standard Model(SM) it is possible to account for the 
total branching ratio when the long-distance contributions are properly included, but all the 



three events are around 214MeV cannot be explained |12| |. If no new evidence support 
the SM explanations in the future experiments with more events, it is most likely to interpret 
the three events with the existence of a new particle, X, beyond the SM. However, a new 
particle explanation for the HyperCP events seems too radical because there are not earlier 
experiments which observe such a 214MeV new particle. If this new light particle does 
indeed exists and contributes to the hyperon decay, it may also contribute to the kaon and 
B-meson decays. So, the fact that lots of experiments at this low-energy region did not 
observe the new 214MeV particle means that strong constraints has been imposed on the 



new particle explanation 
The authors of Ref. 



"or the HyperCP events. 

Sj proposed a argument to explain the HyperCP events for the 
hyperon decay with the new particle X without contradicting with the constraints on the X 
from the low-energy experiments. As shown in Ref. [8[], in addition to the flavor-changing 
two-quark contributions, there are also four-quark contributions arising from the combined 
effects of the usual SM \AS\ = 1 operators and the flavor- conserving couplings of X, which 
are comparable with the two-quark ones and cancel sufficiently to lead to suppressed rare 
kaon decays rates while combining the above two kinds of contributions yields S —>■ pfi^fi~ 
rates within the required bounds. 

Based on the analysis in Ref. [8|], the authors of Ref. p] pointed out that a light pseu- 
doscalar Higgs particle in the next-to-minimal supersymmetric standard model (NMSSM) 
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15| can be identified with X. In fact, the mass of the light pseudoscalar Higgs par- 



H 

tide in the NMSSM can be as small as 214MeV in the large-tan/? limit. Under some 
assumptions, it has been shown in Ref. Q] that there are regions in the parameter space 
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where can satisfy the following constraints: 

B{K^ vr^A?) < 8.7 x 10"^ 
B{Ks ttM?) < 1.8 X 10-^ 

B{B < 8.0 X 10~^ (1) 



which are obtained in Ref. [8| from the measurements of the kaon and B-meson decays 
3, 0, Q, Q, |20| . and simultaneously explain the HyperCP events. 

However, the author of Ref. [3] only considered contributions from the SUSY charged 
current, i.e., contributions arising from the exchanges of chargino and squark, and do not 
include contributions arising from SUSY- flavor-changing neutral currents(FCNC). It is well 
known that the SUSY-FCNC couplings can yield important (and, sometimes, even dominate) 
contributions to low-energy flavor physics, so further investigation on the possibility of the 
SUSY-FCNC mediating HyperCP events is needed. In this paper, we show that the SUSY- 
FCNC effects also can explain the HyperCP events and satisfy all the constraints in Eq.([T]). 



We adopt the mass insertion method 
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24l ] to parameterize the flavor-changing 
effects and calculate SUSY-FCNC contributions to branching ratio of S ^ pA^ and rare kaon 
decays. This method introduce the super-CKM basis for the quark and squark states. The 
couplings of quarks and squarks to the neutral gauginos are flavor diagonal, while the flavor- 
changing SUSY effects are exhibited in the off-diagonal terms of the squark mass matrix 
denoted by {A^j)jj, where I, J = L, R and i,j = 1,2,3 indicate chiral and flavor indices 
respectively, and q = u,d denote the type of quark. The squark propagator is then expanded 
as a series of = (Afj)/j/m^, where m is an average squark mass. Using the mass 

insertion method, we can perform calculations of the SUSY-FCNC contributions to S — > pA^ 
and rare kaon decays. Since the relevant does not involve in the B- mesons decay, we do 

not consider the constraints from B- mesons decay. It is well known 23|, |2J, |25|, |26|, |27|, |28|. |29| 
that the parameters {Sf2)ij used in our calculations also yield important contributions to the 
— mixing, however, our calculations will show that the measurements of the — Ks 
mass difference and the indirect CP violation observable eK do not lead to more stringent 
constraints than ones from kaon decays. 

We organize our paper as follows. In Sec. [Ill we give a brief summary of the NMSSM. 
In Sec. mil we calculate the two-quark flavor- changing contributions to the E and kaon 
decays arising from the SUSY-FCNC effects mediated by neutralino and gluino. In Sec. [IV] 



we combine our two-quark contributions with four-quark contributions in Ref. 8j to give a 
numerical results and discussion. Feynman rules and analytical expressions for the four- 
quark contributions are collected in the Appendix A and B, respectively. 



II. NMSSM 

In order to make our paper self-contained, we start with a brief description of the NMSSM 
and the relevant couplings considered in our paper. The superpotential of the NMSSM is 
given byQ, Q, Q 

W = QYuHuU + QYdHdD + LYMdE + XH^HuN - (2) 

where Hu and Hd are the SU(2) doublet with the hypercharge 1/2 and -1/2 and are re- 
sponsible for the up- and down-type quark mass, respectively. The ratio of the vacuum 
expectation values(VEVs) of the and Hd is defined as tan (3, which are just like those 
in the the minimal supersymmetric standard model(MSSM). Compared with the MSSM, 
there is one more gauge-singlet Higgs Field with the hypercharge and VEV x in the 
NMSSM. After breaking of the supersymmetry, there are seven physical Higgs bosons in 
the NMSSM, including two charged Higgs bosons, three neutral scalar and two pseudoscalar 
Higgs bosons. 



The Higgs potential of NMSSM is|30| 



VHiggs = Koft + Vf + Vd, (3) 

where 

Koft = mljHdl' + mljH^\^ - {XA^HAN -hA,N^ + R.C.), 
Vf = \X\WHd\^ + \H^\^)\N\^ + \XHdHu-kN^\\ 

Vd = ^^(I^.P - \Hu\') + ^^\HiHd\\ (4) 

The above Higgs potential has a global f/(l)_R symmetry in the limit of vanishing parameters 
Ak, Ax 0[31]. If the global U(1)r symmetry is broken slightly, the lighter pseudoscalar 
A\ has a natural small mass: 

,4.=3fa.4. + o(J^), (5) 



In the large tan/3 limit, m^o can be as low as ~ lOOMevjsol [sil- 

The Lagrangian describing the couplings of to the up- and down-type quarks and to 
the leptons are given by Q 

iA° 

^Aqq = -{lumuU-i^u + Idrridd-f^d) — - (6) 

V 



and 



Cm = ^75^1, (7) 



respectively, where 

Id = -9i = ( .T^'" ) ' = — (8) 
a/2x \Ax + kx J tan^ /3 

Note that can be neglected in the large-tan /5 limit. The four-quark contributions can be 

deduced from the interactions in Eq. ([6]) combined with the operators due to W exchange 



between quarks j?]. 

It has been shown in Ref. 7|] that an A\ of mass 214.3MeV decay dominantly to muon- 
antimuon pair, and B{A\ Ai+yU^) ~ 1 can be assumed. In addition, the constraint imposed 
by the muon anomalous magnetic moment is given by 

\9i\ ^ 1.2 . (9) 

Moreover, the neutralino sector of the NMSSM is different from that in the MSSM. There 
are five neutralinos in the NMSSM, and the Lagrangian for the mass term of the neutralinos 
can be written asfisl] 

C^^,=-\{i^'fYi,' + H.c., (10) 



where Y is the symmetric neutralino mixing matrix ,its expression can be found in Ref. |15 |. 
The masses of physical neutralinos can be obtained by diagonalizing Y by a unitary 5x5 
matrix N 

m.^o6ij = N*jr.mnNjn- (11) 

III. THE SUSY-FCNC EFFECTS 

The full hadronic amplitudes for the kaon decays in Eq. ([T]) and S pA\ all can be 
written in the following form 

A^fuU = M2q + Miq , (12) 
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FIG. 1: Feynman diagrams for \^S\ = 1 transitions, with I, J = L,R, and i,j = 1...5 

where M.2q arises from the SUSY-FCNC interactions at the quark level, and M-^q arises from 
the four-quark interactions. Using the method shown in Ref. js], we calculate the four-quark 
contributions, and their expressions and numerical results are given in Appendix B. In this 
section, we mainly concentrate on the hadronic amplitudes M.2q from two-quark contribution 
and give our analytical results. In general, there are two kinds of FCNC contributions arising 
from neutralino and gluino exchange, respectively. The relevant Feynman diagrams for the 
s ^ d transitions are shown in Fig. [H Note that there are five kinds of neutralinos in the 
loops in the NMSSM. Calculating these Feynman diagrams, one can obtain the two-quark 
FCNC Lagrangian for s dA\ 

C-Asd = iCid^—^sA^^ + iCfiJ^i^sA? + H.c, (13) 
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with 

CLiR)-Clln) + ClRy (14) 

where Cf(^) and C'f (^) denote contributions from neutrahno and gluino exchange, respec- 
tively. They are given by 

5 

+ i^i2)LR [lljmdmsL2iRlj - I-mdmsL2iRli - lfmdms{LuRli + -t'2i^2i)] 
HK2)rr [ll.msLuRlj - ifmMiRu + ifm^imlL^iRl, - Lr^R*,)] } , 

5 

+ {St2)LR [liR2iLlj + I!R2^LI^ - lt{mlRuL\, + m^i?2«L;j] 
+{^12) RL [lijUidmsRuLlj + ifmamsRiiLli + lfmdms{R2iL*2i + -RijLtJ] 
+ {St2)RR [lijrridRiiLlj + I^mdRuLl, + lfmd{R2iLli - mlRui;,)] } , 
CI = ^C^rn,[2VAddCi{y-M2)RL + V^^^D,{y-^^^^^ 

d 

C'r = ^C^rn,[2VAddC^{yM2)LR-V^^^^D,^^^^^^ (15) 

d 

with 

VAdd = 7^ nU[l, 

Zmw cos p 



cw 



2i 



V2mwSw cos P ' 

^ _ g^^^i,! - 3ctyiVi,2 ^^^^ 

II^ = D2{y^,y,)R^:^ + ^^D2{y^,yi)R'(;,, 

4 = i^2(2/.,%)i?5 + ^^^g^/^2(l/.,%)i?5, 

= '2VAddm^fi^{yi), 
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Uq 9! Ui ^ i...O, 



(18) 
(19) 



where is used to diagonalize the pseudoscalar Higgs mass matrices, the matrix A^^-, i, j = 
1 ■ ■ ■ 5 is the 5x5 unitary matrix defined in SecUll cw = cos6'vi/, where 9w is the Weinberg 
angle as usual, and 



W'R 



-[iU[,cosP + U[,smP) 

P „;„ o ttP 



i4, 



cw 



-V2kU[,{N,,N*+N,,N:,) 



W'L* 



(20) 



And the loop functions Ci{x), -Di(x) and D2{x,y) are defined as 

X — 1 — X log(a;) 



Ciix) 
Di{x) 



{1-xy 

1 — + 2xlog(x) 
2(1 ' 
1 

{l-x){l-y) {l-xy{x-y) (l-?/)2(y-x)' 



xlog(x) 



(21) 



The quark level effective Lagrangian in Eg. (1131) can be mapped onto the chiral Lagrangian 
at the leading order ^ 

Ca = bD{B{hA,B}) + bF{B[hA,B]) + bo{hA){BB) + ^f^Bo{hA) +li.c., (22) 



where 



/ 1 vo I 1 



B 



p 



\ 



7:0 



represents the baryon fields, /^r = 92.4MeV is the pion decay constant, (■ ■ ■ ) = Tr(- ■ ■ ) in 
fiavor-^f/(3). And 

hA = -t{Cneh^^ + CL^hOA'i (23) 
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where 



'^O 0^ 



h = Tq + iT-j 



1 




is used to specify the s ^ d transition, 



(24) 



and 



TT 



1 

7! 



\ 





7r+ 


K+ ^ 








K- 







is the pion octet. 

The two-quark amphtude A42q can be deduced from £^ plus the usual chiral Lagrangian 
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34 



35| 



Cs for the strong interactions of hadrons, which is expressed as 

= t{BYD,B) -mo{BB) + DiBYldA^B}) + F{BYl5[A„B]) 
+bD{B{M+, B}) + hF{B[M+, B]) + 6o(M+) {BB) 

with 

D^B = d'^B + [V'',B], 
where M = diag(m, m, rris) is the quark mass matrix in the m„ = = m limit. 



(25) 



Using the mass relations ms — nip = 2(6/) — 6i?)(ms — m), m 



m: 



Bo{ms — rh) and 



m|; = Bo{ms + rh) , the amplitudes for the different decay modes can be written as 



(26) 



+i{D - F) 



Cl-Cr\ 5o(ms + rrip) ^ 



where Bo = 2031MeV|4| 



(27) 



IV. NUMERICAL RESULTS AND DISCUSSION 



In this section we present our numerical results. The Higgs sector of NMSSM is describe 
by the six independent parameters 

A, k, Ax, Ak, tan/3, /i, 

where fi = —Ax. For convenience, we will take m^o and the coupling of down-type quarks to 
A^, Id instead of k and A. We follow Ref. [?] to set la = 0.35, —Ax = ISOGeV and tan (3 = 30. 
y4fc and Ax are set as 0.001 and 0.002, respectively. The mass of A^ is set as 214.3MeV to 
satisfy the HyperCP data. We set the mass of gluino and average down type squark mass as 
200GeV and 350GeV, respectively. With our input, the mass of the neutralinos are around 
100 ~ SOOGeV. Numerically, contributions from the exchange of the gluino and squarks are 
lager than contributions from the exchange of the neutralinos and squarks, i.e., C^°^ are 
larger than in most regions of parameter space. This is due to the effects of as- 

Our numerical results are shown in FigsIS] - [71 We first assume that are real. 

The allowed regions in parameter space are shown in Figs|3] - [5l where the grey areas are 
the allowed regions for A^ to explain the HyperCP events. When the constraints obtained 
from the kaon decays are considered, the allowed parameter space are greatly reduced to 
the dark regions. From FigslH] - (Sj the constraints on the combinations of the parameters 
can be obtained as following: 



{^12} LL{RR} iSf2)LR{RL) 


< 


3.9 X 10' 


-12 


(without kaon bounds), 






3.7 X 10' 


-12 


(with kaon bounds). 


{^12) Ll{8^2) RR 


< 


3.9 X 10" 


-9 


(without kaon bounds). 






1.8 X 10' 


-14 


(with kaon bounds). 


(^12) LRi^u) RL 


< 


2.7 X 10' 


-16 


(without kaon bounds). 






0.9 X 10' 


-17 


(with kaon bounds). 



It has been widely studied in the literature that the SUSY-FCNC effects has great impact 
on the — mixing if the relevant (5i2)/j are complex. So we further investigate the 
possible constraints on {Sf2)iJ from the Kl — Ks mass difference ArriK and indirect CP 
violation parameter ex- However, the constraints shown in Eq. (l28|) are roughly several orders 
smaller than those given in the literatures involving the SUSY-FCNC mediated — 
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mixing |23|, ^ , l25|, l26|, l27|, l28|, l29| , where (^fs)^^ and {5f2)LR{RL) are around 0(10"^ ~ IQ-^) 
and (9(10^^ ~ 10~^), respectively. This fact indicates that the constraints from the — K^ 
mixing may be automatically satisfied once the constraints from E pA\ and the rare kaon 
decays in Eq.([T]) are satisfied. Our numerical results do indeed confirm that the constraints 
from the — mixing do not lead to more stringent constraints than those ones given 
from the kaon decays in Eq.([T]) . 

Figs|6] and [7] show the constraints on the complex {5f2)ij from E — * pA\ and rare kaon 
decays in Eq.([T]). And the corresponding constraints on the combinations of parameters are 
given by 

R'e(5f2)Li?(i?L)Im(5f2)LiJ(i?L) < 2.7 X 10~^^ (without kaon bounds), 

1.5 X 10~^^ (with kaon bounds), 

^^{^12) LL{RR)^'^{^i2)LL{RR) < 1-8 X 10"^^ (without kaou bounds), 

1.8 X 10~^^ (with kaon bounds). (29) 

From FigsIH] and [71 it can be seen that the grey areas are the allowed regions of the SUSY- 
FCNC parameters for A\ to explain the HyperCP events and the grey regions are greatly 
reduced to the dark ones when the constraints from rare kaon decays are considered. Even 
so, there are still regions in the SUSY-FCNC parameter space where A\ in the NMSSM can 
be used to explain the HyperCP events without contradicting with the constraints from the 
rare kaon decays and the — mixing. 

In conclusion, we have calculated the two-quark contributions to the decay S"*" p^~^fi~ 
arising from the transition s dA^ via the SUSY-FCNC couphngs. Combining the two- 
quark contributions with the four-quark contributions, we show that there are regions in 
the SUSY-FCNC parameter space where the A^ in the NMSSM can be identified with a 
new particle of mass 214.3MeV, X, which can be used to explain the HyperCP events, while 
satisfying all the constraints from the measurements of the rare kaon decays. And once the 
constraints from the kaon decays are satisfied, the constraints from the — mixing are 
automatically satisfied. 
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Appendix A 



We give the Feynman rules used in our calculations in FigJJl where V^jj, VAdd, -^(-R)i{2)i 



and Riif^^ are defined in Eq. f|T6|l . Eq. f|T71) and Eq. fl20l) . respectively. And 

Sij = 




Appendix B 



We collect expressions of the four-quark amplitudes for the different decay modes in this 
appendix. A detailed description can be found in Ref. {s], we cite their results here. As 
pointed in Ref. [7|, the couplings of to the up- type quarks tends to zero at the limit of 
large tan/3, so we neglect terms that are proportional to /„. 

The four-quark contributions for the kaon decays are as follow: 



se\ 



K 



ml)se 



+ ^8(7 



6(m2 



+ 



{2m\ + ml - 3m\o)s0 



{Am\ - ?>ml)se - V2{2m 



K 



ml)ce 




(30) 



{Ami - ^ml)ce + V2{2ml 



ml)se 



6V2( 



+ 



m 



a; 



mt 



+A/8(m|- - ml)ce 



(4m^ - ?,ml)s0 - y2(2m2 



6\/2(m^o - "^5') 



(2m|- + ml - 3m^o)s6i 
(31) 



K 



ml)ce 



12 



dj dj 





\ 
I 



I 



d^ x° 



I rff I dl 



i^9.,{5niPL - SilPM -ie[{Lu5ii^ + rriaLMPL 

+ (nidRiiSiR + R2i5il)Pr 



d d dj di 

'a? 'a? 

I I 
I I 

FIG. 2: Feynman rules used in our paper. 

where 73 = —7.8 x 10~®, sg and cg are short for sinO and cos^, 6 = —19.7°. And, the 
four-quark contributions to the S pA^ process are 

A^4,(S+ ^ M?) = ^P{ApAl - 5p^o75)E+, (32) 

with 

^pAO — I'd 77— \ 2 2" 2 2 

^ (4m^ - ^ml)sl - V2{2m\ - ml)ceSe | ^^^^ 
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and 

D , U Bp^o { ml {Am\ - 3ml)cl + \/2(2m^ - ml)c0Se 

^pA° — I'd — \ 2 2 2 2 

where A^^o = -3.25 x 10"^, Bj^o = 26.67 x 10"'^. 
Numerically, the above amplitudes are 

M4g{^ pA?) = ip(-6.96 X IQ-^ld— - (5.71 x 10-^)^—^5)^+, 



m 



K 



Miq{K+ 7r+^°) = -il.08 X IQ-'^ld 

u 

2 

TYl 

MAq{K^ ^ ttA?) = ^1.12 X (35) 
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FIG. 3: The allowed values of (5i2)iR = {^12) RL as a function of {^12)11 = {^i2)rr- The grey 
area is the regions where can explain the HyperCP events, when the constraints from rare kaon 
decays are considered, the allowed regions are greatly reduced to the dark ones. 




FIG. 4: The allowed values of {5i2)ll as a function of {Sf2)RR, The grey area is the regions where 
Al can explain the HyperCP events, when the constraints from rare kaon decays are considered, 
the allowed regions are greatly reduced to the dark ones. 
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FIG. 5: The allowed values of (5i2)i?L a function of {§12) lr. The grey area is the regions where 
Al can explain the HyperCP events, when the constraints from rare kaon decays are considered, 
the allowed regions are greatly reduced to the dark ones. 
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FIG. 6: The allowed values oi ImliSf^) lr = iSf2)RL} as a function oi JieliSf^) lr = {Sf2)RL}- The 
grey regions denotes the survival regions for explaining the HyperCP events alone; In the dark 
regions, the can explain the HyperCP events and simultaneously satisfy the bounds originate 
kaon decays and the — mixing. 
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FIG. 7: The allowed values of Im{((5f2)LL = {Sf2)RR} as a function of Re{{5f2)LL = {Sf2)RR}- The 
grey regions denotes the survival regions for explaining the HyperCP events alone; In the dark 
regions, the can explain the HyperCP events and simultaneously satisfy the bounds originate 
kaon decays and the K'^ — mixing. 
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